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ABgERfcgP 


■-/ik Ifonte-Garlo siiaulation tedmique suggested by 
Shinozuka is used to investigate tie response characteristics 
of a few simple lin^r and npn-linear structures, under 
random vibration environment. Stationary response characte- 
ristics are computed for both linear and nonlinear systems. 
Jov the linear system laie transient response is also computed 

Results of simulation are used to calculate the 
first-passage probability density, the ^pected amount of 
time spent above a prescribed threshold, the number of 
crossings of a level per unit time. A linear single degree 

of freedom systan is considered for solving the above 

/ 

laeii'tiojied problems* v 

The random walk model developed by Toland is exten- 
ded to a two degree of freedom systoa, and the computation 
problems Involved in solving such a mod^ are discussed. 
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GEAPIEH I 

mmmmm 

1.1 In the past two decades there has slgniJtlcant 

devilopmen-|f In the theory of Stochastic Processes. This has 
na'de it possible to predict idie djneanic stame-fajial response 
nhder ^ndom loadli^ conditions. It has been studied erfe enslvely 
wildti Id^ly idealised load structure systea®.. Howerer, the real 
life loadix^ conditions and structiiral systene differ signifi- 
cantly from the assumed deteaninistic models. lioadii^ conditions 
and materiaG. properties in almost all real life structures have 
a large measure of imcertainity requiring a stochastic approach 
for areallstic treatment. Hence a probabilistic, approach, pro- 
vide a more realistic approach to assess the overall structural 
performance. 

The probabilistic approach too has its limitations, 
loading conditions, 'many a time, are nonr-stationary in nature. 
The well dev^oped principles of stationary sto<hastic processes 
cannot be applied to ir»nstationaiy stochastic processes. Also 
it is limited to the linear structural models. Attempts have 
been ma^le to solve some of the problaa©, which resulted in the 
equivalent linearization techniques and perturbation techniques. 
However, these are restricted to structures wilh weak non- 
linearity. Some of the problems such as the first passage 
probability e^nd fatigue damage of a structure, which form the 
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more important failtjre criteria^ are yet unsolved* Monte Cfetrlo 
technique have onerged out as powerful tools and have heen 
used successfully to solve the above mentioned parcblems. 
lEhese are essentially numerical techniques. Some of the problons 
which are not amenable to cozwentional mathematical amlysis 
are* a) analysis of panel vibrations of aircraft and sutsBarines 
inAiced by boundary layer turbulence* b) analysis of off shore 
structures taking the wave, structure interaction into consi- 
deration, c) analysis of vehicle and aircraft vibration caused 
by the roughness of the ground etc. 

The Monte CJarlo techniques basically consists of simu- 
lating the loading conditions, which have fixed date and solving 
the resulting equation of motion using numerical integration. 

The artificial data have the same statistical properties ag 
that Qf parent one* Por a stationary stochastic process a single 
sample response function is enou^ to compute the statistics 
such as mean, variance etc,, as the principle of ergodicity 
is valid. The temporal averages can substitute ensemble averages. 
Wherever the stochastic process of interest is nonstationary, 
an ensemble of sample fmctions can be collected and the averages 
can be comput ed across the ^semble. Thus by using artificial 
data for the loading conditions statistics such as peak response 
distribution, number of times a critical value is attained, 
transient bdiaviour of the structures, can be conveni^itly 
con5)uted following the Monte CSarlo approach, .The simulation 

' w 
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techniques have the advantage over the convent ioml methods 
in that these can he applied to any type of load structure 
syst ems . 

Ohe limitation of Monte Carlo approaches is that 
often they require a consideiahle amoxmt of comtaiter time* 

Analysis of noivstationary vibrations require large sample 
size* particularly for the structures with low damping, it 
requires a very large number of sample functions. Por Ihe 
existing computer system, IM 7044, the time consumed is of 
the order of a few hours. However, the advent of third gene- 
ration computers should overcome this problm. 

1.2. Literatur e Survey. 

1,2.1 3Jhe basic principle of stochastic process theory 
and its application to the response ai^lysis of structural 
systans have been discussed ertenslvely by bin Q] and 0randall[2] 
In their worfcS-Bjth tine and frequency dDoaln approaches are 
elaborated to deal with the response charaotorlstios. But these 
methods beoooe extraaely oombursome and involved when dealing 
with the response characteristics of oomplloated systems 
especially subjected to noh- stationary excitations. Honte-Carlo 
simulation teohniq,ues have been proved to be extranely useful 
in dealing with such cases 

Filtering techniques were being used hy mny authors 

before Borgman [s] -published an efficient means pf simulating 
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a sample function of a stochastic process. He ccmsidered 
the ocean surface elevation as a random process ai^ rep^csa^ited 
it as a sum of a series of cosine functions with random phase 
angles and amplitudes consistent in the energy in the target 
spectral density at that frequency. 

Later Shinozialca [5,43 proceeding in the similar lines 
has developed a more generalised approach to simlation* which 
can he applied to multi-variate and multi dimensional processes.^ 
Also it has been shown that among all possible series expansions 
of a random process, Shinozuka*s approach gives the best 

approximation, 

A number of such diverse problems lihe, the nonlin^r 
panel response from a turbulent botindary layer [Ylr J^i’^linear 
dynamic analysis of off shore structures K, response of the 
structures to wind loads have been solved receitly following 

the Monte-Carlo simulation technique of Siinozhka. 

1. 2,2 Transient Response Analysis 

The transient response of a system under landom 
excitation has a h^-story extending over the past 40 years. 
Dhlenbeck and Ornstein gave the first solution to the problem 
in 1950. Later Hil^beck and Wang ^5] have solved the problem 
using the fokker-Planck equation for «ie system* Sibsequently, 
several authors have considered the response of dynamic^^ 
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systems confi.tding tiiaaselves to the statix>mry aspects of 
the motion, Causey and Stun^tf jji’J have analysed first, the 
transient motion of a sin^e degree of freedom syat®i, subjected 
to a stationary random input process, having an arbitrary 
power spectrum. In the same reference an application to the 
response aialysis of a structure subjected to strong motion 
earthq.uafce is given. Since that time, a significant amount 
of work has been contributed by sevei®.! authors, mainly in 
the area of earthquake engineering, dealing with the non- 
stationary transient, response analysis of the structures 
subjected to earthquake accelerations. However, the simulation 
techniques have dominated, whence « number of sample response 
functions were simulated artifically and the response charac~ 
teristics were computed from the ensemble averages. 

Recently, Chakrovorty [j^ has developed an analytical 
approach following the evolutionary' spectrum concept of 
Priestly 151} , for the non-stationary response analysis of 
linear structures subjected to the white noise excitation. 

In. the same reference, an example is given where the trapsi^t 
motion of an equipment in a building which subjected to earth- 
quake acceleration is studied, The support acceleration of 
the equipment becomes non-stationary because of its transient 
nature. 
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1*2*3 Pirst Passage £roM«wi p 

Urst exciirsion failinre pxobaTjllity and the fatigue 
damage are two JLmportant modes of failure in the area of landom 
Tihratio3as« !Ehe former relates to the time at which for the 
first time the response of a structure crosses a specified 
threshold and the later to the cumulatiTe damage during 
vibration* So far Idle exact solutiop?' even for the simplest case 
that of a narrow band random process such as the output 
of a li^tly damped linear oscillator, has not be established* 
However^ based on the numerical results certain analytical 
approaches have been ^i^ested which could establish ihe 
bounds for such probabilities 4 pollowing is a brief survey 
of the previous work in, this area* 

Some of the first passage problems for a li^tly damped 
linear oscillator excited by a white noise have been solved 
by Crandall, Ghandiramani and Cook {j!^ in 1966, They described 
two naimerical approaches to the problem. In one ihe probability 
statistics of safe intervals were made use of to solve the 
problan and in the other a markov vector approach was used* 

Two types of initial conditions, the zero start and the stationary 
start and three different lypes of ^fe domains the single 
hajnrier, the double iKi-rrier and ihe envelope type of barrier 
were considered* 

Crandall [143 described variety of analytical 
approximations to the first passage probability density for 
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linear oscillator and con^jared with, results obtained by 
simulation and numerical work* He derived expressions for the 
first passage probability density with iiie following assumptionsi 

1 ) Assumption of independent crossings 

2) Assmptions of independent peaks 

3) Assumption of independent envelope crossings 

4) Assumption of independent envelope peaks 

5) Two state Markove Process assumption, 

Ariaratnam and Hennon [t?] have solved the first 
passage problen associated with snap - throu^ of a shallow 
cylindrical ^ ell subjected to a wide band stationary stochastic 
loading. They followed the numerical approaches given in 
Ref. O 33 . 

Recently Toland, Ife-ng and Hsu [l 6, 20 have proposed 
a random walk model for the linear oscillator and solved the 
first passage problem. The details of the model are discussed 
in the CSiapter Y. 

1.5 The Present; Thesis Work 

In the present thesis work the response charactacis- 
tics of a single degree of freedom linear system to band 
limited white noise excitation are conqputed to illustrate 
the Monte Carlo simulation techniques. As an application it 
is extended to the response analysis of an aircraft— landing 
gear assembly to taxiing induced vibrations due to runway 
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undu tat ions. Briefly the contents of the thesis are the 
followingi 

A simtilation procedure to generate sample functions 
of any statioisiry gaussian randcm process, with prescribed 
statistical properties is presented in Cle-pter II. Chapter III 
contains an analytical trealment, which makes the lise of 
evolutionary spectrum concept, for the analysis of non^-stationary, 
transioit response of the linear systan, followed by a numerical 
approach. Also a Monte-Oarlo approach to compute the statioiE.ry 
response characteristics of a sin^e degree of freedom systaa, 
an extension to a nonlinear system of Buffing’s type and the 
stationary response ai^lysis of an aircraft -landing gear assembly 
dxiring taxiing and take off are given. In Ihe Chapter W the 
first passage probability density and other related stetistics 
such as the expected amount of time spent above a prescribed 
level, the expected number of crossings of a level, mean time 
to failure are discussed. A random walk model for 1he ar^lysis 
of first excursion failure probabilities is presented in 
Chapter T, 



chapter II - 

■fc., 

SlMOLkHOH 

II. 1 ^st structural systems operate in random envirorment. 
lo predict the dynamic structural performance, therefore, either 
a sufficient number of sample functions of the mvironment should 
be available or such sample functions may be simulated based 
on available stochastic properties. This chapter presents a 
brief discussion of the simulation procedure employed for the 
present thesis work, 

II. 2, Two basic approaches are available to simulate a random 

process: 1)» ®3.e required random process is obtained as the 
output of a linear syston having an appropriate frequency 
response function and is modulated by an envelope function. 

G3iese filtering techniques dominated for quite a long time, 
though they are limited to a one-irariate, one dimensional 
random process. 2) The second approach is based on super- 
position of a series of sinusoidal pulse shapes each having 
a random phase and an amplitude consistent with the target 
power spectral density. It is Borgaan [Sj who first simulated 
the ocean surface elevation following the above approach. This 
has an advantage over the previous one in that it can be 
extended to multivariate, multi dimensional random processes 

also, hater Shinozuka . [3j, £3, dev^oped the wave super- 
position method for simulating both stationary and non-statiomry 

Random Processes. 
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Shijwzulca’s method has been followed in the present 
work. I!he governing equations are as follows s 

k homogeneous, Gaussian, one variate and one dimen- 
sional random process, 3^t), with zero mean, and spectral 
density function^ So (w) can be expressed in the form of 
the sum of the cosine functions by 

r 

x^(t)=T2 S oos (Wj^:t - <)j,) (II.1.1) 


where are random variables distributed uniformly between 


0 and 271, and are iindependoit',- 

K - fo. ) AW J 


(II. 1.2) 


(k-l/2)Aw 


(II. 1.3) 


Gr^(w ) is the one sided power spectral density and 
is related to the mathematical two sidei-f.S.D. by 


G^(w ) = 2 Sq(w ) 

A sample function i(t) of the random process x(t) can be 

rated digitally on a computer, by replacing the random 
variables ^ by their realised valu^. 

U 

i(t) ss f2 S Aj. cos (II.1.4) 

k»1 

The normality of the simulated process is ensured by 
virtue of eenfcral limit theory |^lj . It has been shown in 
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Ref . [4j liiat tlie temporal and ensemble averages coincide, 
at least upto second moment, in limit to the target values* 

This establishes the ergodicity* This makes the method directly 
applicable to a time doc^in analysis where the ensmble 
av^ages can be evaluated in the form of tanporal averages. 

II. 2. Tests for Stationarity ; 

In this section two ron-parametric statistical tests 
are described by means of which the stationarity of the 
simulated random process can be tested* 

Once a sample function of a stationary random process 
has been generated, it is of censid^able interest to know 
how far the simulated sample deviates from stationarity and at 
what level of significance the data can be trusted. Stationarity 
implies that all statistical properties of the stochastic 
process of interest are invariant with respect to tim^ tran- 
sactions, ^der follcwing assumptions, practical tests for 
testing stationarity can be developed. 

1) lor a stationary stochastic process, if we 
compute the statistical propiKcties such as means and m^n 
sq.uare values, by time av^aging ^'|5Arer a sequence of ;^rt 
time intervals of a sample record, the values would not vary 
too significantly from each olSuKP. 

2) It is enou^ to test for w^k statiomrity 
verifying strong stationarity. 
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5) IHie length of the sample record is very long 
con?ared to the random fluctmtions of the data time history, 

4) If the mean square value is stationary then the 
auto correlation function of the data is also stationary. 

The sample record is divided into a number of 
intervals and the mean and mean square values are computed 
to form sequences as follows: 

^2^ ^3 • • • ^n 
-2 -2 -2 -2 

The above sequences are tested for underlying 
trends or variations other than due to expected sampling 
variations. The tests are described below. 

Run Test { 2 ^ 1 - k run is defined as a s^uence 
of identical observations that are followed or preceded by a 
different observation, or no observation at all. let the 
sequence of measured values x^, i=1, 2..*Nhas a mean x. 
Suppose we compare each x^^ with x and form a sequence of plus 
and minus observations, plus indicating that- x^ ^ x and minus, 
< X, we mi^t get a sequence as follows: 

+ + 4-1 ■ 4 -' . , . . + + 

t 2 3 4 N 
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From the sequence the number of runs is counted as indicated. 
The number of runs occur in a sequence gives an indication 
whether the values are independent or biased. The distri- 
bution function of rCk) is presented in Table 4.11, page 170, 
of Ref.[2|]. 5?be hypothesis that the simulated sample function 
is stationary can be tested by observing the number of runs 
in the sequence of means, and mean square values at a desired 

level of significance. 


Trend Test ?- 

It is similar in application to run test. Suppose 

^ ■? — 1 N The number of times that 

we have a sequence x^, i - 1, 

X > x^ for i<j is counted. Each such inequality is called 

arrangement. The total number of reverse arrangements 

is dsBC'fcBd by A* 

k is defined D^tbematically as follows. Eor the 


set of observations x.|, X 2 ...x^, define 

- 1 > ^3 i<3 

0 otherwise 


then 


and 


N 


A. = I 
1 3=i4-1 


H-1 

it *= E 

i»1 ^ 


(II.2.1) 


After counting the number of reverse arrangements 
la a sequence, tbe hypntbesis can be tested at any dnaired 
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level of significance, by cbecking ?3iether it is falling 
within th.e corresponding limits specified in the standard 
tables liable 4*12, page 171» Ref.{22j]. 

Power Spectral Itensity Calculations 

In this section, a method of estimating the power 
spectral density, and choice of various parameters involved 
are discussed. 

Def. Power Spectral Density of a stochastic process 
is defined as the Ponrier transform of the auto correlation 
function as follows. 

G-^(f) = 2 / i^Cr) e dr {II.5.1) 

— ..oo 

where ® x(ti- r)3 

6^(f) is one sided power spectral density and is related to 
the mathsiatical two— sided power spectral density as follows. 

s^(f) - S^{-f) = a^(f)/2 (II.3.2) 

In order to estimate the power spectral density from 
the simulated sample function of the random process of 
interest, Eg.* II.3.1 shoiild be discretized. It can be easily 
shown that [2^ 

= 2h IS, + 2 2^ Sj c«s (rok/m) +(-1)*^ . 

(I1.J.3) 
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The parameters h, m and fc should be chosen appropriately 
to get a good estimate of P*S.D.. function. It can be j^own 
that the normalized standard error involved in the estimation 
of P. S,D. from Eq.n, 11.5*3 is given by 

» m/N (11,3.4) 

where N is the sample size. 

The maximum lag number m, is also related to the resolution 
band width Be of the P.S.B. estimate by 

mh * 1/Be Cll.3.5) 

Thus a large m, thou^ increases the resolution of the 
P. S.D, estin^te, will increase the uncerfcainity of the estimate. 
So a compromise should be arrived in selecting the value of m 
such that sufficioit resolution with minimum statistical 
uncertainity is ensured. Usually m is chosen less than one 
tenth of the sample size N. 


gives a raw estimation of the P, S.D, function at dis- 
frequency intervals (^^), k = 1, 2, ...m, 

fc is cut off frequency and is given by fe = 

h is sampling interval 

is estimate of correlation function at llag r 
m is mxiDKBa lag number 
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Tile auto correlation function is estimated by 


^ “ sir * 


'OkT 


(II. 3 .6) 


ss 0^ "ty 2f #iii 


As the value of h is usually fixed before hand, the value 
of «utc-<iff. frequency, frequency beyond vdiich the P.S.D. 
has no significant value, should be set at least eqml to 
l/ai, which is called the Nyxjuist frequency, 

fc > l/2h (11,3.7) 

G (f) will be only a raw estin^te and it sboiiLd be 
smoothed to get a better estin&te of the true P.S.b, Ihis can 
be done by multiplying the auto correlation function with a 
‘lag Window*, which tantanK) tints to weiring the correlation 
function non-uniformly. The final formula is 


A 


= 'g' (— ) = ax % cos 


(11.5*8) 


where 


1 + cos -- 


r = 0, t, 2, /..m 
. r > m 


(11.5.8(a)) 


for Hanning* s lag Window, 


and, 


3 2 , 
1+6 (|) - 6(rAl) , r « 0, 1. 2, ...m/2 


D - (|a' 


9 ^ 2 


- 2ti» • • • 


m 


* 0 


r > m 


(11*3 .8(b)) 
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for Parzeii*s Iiag Window. 

Any one of the above two lag windows can be used 
satisfactorily, but it is found from the results that for 
power spectral density estimates, Parzen^s lag window gxves 
a smoother and better estimate than Hanning's lag window. 
The effect of varying m is also studied. It is found that 
as m incr^ses the scatter in P.S.P. estimate increases. 
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JJ,/} Begin tfi and Discussion: 

Presented Herein are tHe results of the simulation 
technique desorlhed in the previous paragraphs. Two examples 
are considered, ^or the first one a wide hand random process 
with constant spectral density is chosen. Ihe target power 
epectxal density function is given hy the following equatron. 


GqCw) = 4- luiits 


< w < W2 


(II.4.I) 


where 


w = lower liait for frequency - 0 cps 

^ = cut-off frequency = 

2 

No. of terms considered for the simulation are 100. 

Second example is that of runway rou^ess for which 
the power spectral density fihcticn may he written as 

0 ( 11 . 4 . 2 ) 

fL^ <n<il2 

JX = wave nunber in rad/ft. 

s ,CX)1791 rad/ft. 


vhere 


9 is introduced in the equation II.4.2 to avoid the sto^- 
at the origin and to the power spectral densx.y 

= a constant varying from . 5 x 10"® to 20 x 10 


0 
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Ai = 0 

St 2 ^ l/2h-* 

where h is the step length used for the integration of the 
equation of action. 

sample functions are generated according to the 
Eqi^tion wliere 

Wjj: = (k-1/2) ^ w + w, 

K* 1, 2 . . . 

and - (^2 “ 


are generated uniformlj between '0' and '2 b' naking 
of the built in subroutine EBBH of IM 70U systeh. 

In Table So.II.1 the results of tests for stationarity 
are presented. The nuaber of runs and the number of reverse 
arrangements oleo-rly fall within the respective acceptable 
regions thus establldilng the stationarity of the generated 

sample* 


Results of power spectral density estimates are pre- 
aented in Pigs. II.1 II.2. Table Ho.II.2 gives a com- 

parison of the power spectral density estimates from 
function, using Hanning's lag window and Parsen's lag 

It may be inf ered from Jig. II. 1 tlaat 
window. . - -■ - • J 

■inriow eives a smoother estimate. Prom the 
i^rzea’s lag window gives 



20 


Fig, II, 2 ii cQ-n be seen thst tiie scatter in ttie estinsite 
of jHDwer spectDral density increases with, lag mmber m, 
which follows from Eq. 11,3.2. Eor n equal to 20 the hounds 
are 4+. 283 2 and for h equal to 80 the corresponding limits 

are 4 + .5684, 
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RESPONSE? MiJiYSlS 

Any structural system wh^ subjected to an excitation 
undergoes a period of transient vibration. If the excitation 
is a stationary lendom process, the structural response even- 
tually becomes stationary* However, the tiansient respoime 
is non-s tat ionary, IQiis chapter provides a brief discussion 
of an analytical approach which uses the concept of 

evolutionary spectrum- to find the transient response of a 
simple linear oscillator subjected to white noise excitations 
and a numerical technique based on the ensenble averages, a 

Monte-CJarlo technique for finding the statiomry response 
characteristics of a linear and nonlinear sin^e degree of 
freedom systems and an idealised model of an aircraft ~ 
landing gear assembly* Nmerical approach utiliEHgsthe 
results of simulation, 

2. 3516 characteristics of a linear time invariant system, 

which is shown in Hg *311^1 r c£m be represented by the well known 
differ ^tial equation 

y^t ) < 1 - 2 ^ Hy( t) + w^y(t ) * x(t ) (I H *2, 1 ) 

2 , 

wh©ce V4 * k/m 

^ « c/2mWo 
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Wo is the iiatural frequency , c is the Tiscous damping coefficient, 

^ is the fraction of critical damping, x(t) is the base exci- 
tation in terms of acceleration and y(t) is r Native displacement, 

Ihe auto correlation function of *:x(t) may be written as 

= EpeCt^) (111.2,2) 

Ihen the process is weakly stationary the auto correlation 
function can be written as 


EjCi: ) * / %(w) expliwCtg-tpJ dw (III. 2.5) 

where 

X = tg - t^ and S”(w) - two sided power spectral 
density of input Random Process, 

Ihe response y(t) is non- static nary, as it builds 
up from rest, and its auto correlation function may be written 

RyCt^ftg) = E|y(t^) y(^^ (III. 2.4) 

The response of the oscillator y(t) is given by the 
(jojpmlution of the impulse response function with input exci- 
tation x[tj. Thus, .i 

y(t) - / iSCt) h{t- r) dt (I1I.2.5) 

where h(t ) is the im^se response function and h(t-t ) * 0, . 

:for ;t < 



Substituting ©i, III.2,5 in Eq. III.2.4, 

SyCt^jtg) ~ E|^ J" x( ii(t^*“ 7*) d Z.^! f x( ^2^ 

^ o o 

digj (III. 2.6) 

After simplifying the Eq. 111-2.6, we get the following 
relation for 

K(w»ti) K%,t 2 ) exp|iw(t^-t 2 ) 3 (S/f 

(111,2.7) 

t 

where K(w,t) » / h(u) exp(-iwu) du (111,2.8) 

and E*(w,t) is the contplex conju^te of K(w,t), E(w,t) is 
analogous to the frequency response function and is called 
time d^endent transfer function. 

The impulse response function h(t.) and the frequency 
response function H(w) are givai by 

h(t) « - ^ exp[->^otO ^ t 2 0 

= 0 ; t < 0 (ni..2.9(a) 

HCw) = - (III.S9(b)) 

K(w,t) may be obtained as lo} 

E(w,t) »: H(w) ^l-fcosw^ct sinw^r.t]exp|3rSv'V+i^)‘yJ 

(111,2.10) 
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It is clear from Bi.III-2.10 that K(w,t) -* H(w) as 
t-^. The time dependent mean square response is given hy 
putting t-j = t^ in the Eq« III-2«6. 

S[7^(t)l *=<^y(t) / Sjj. (w) K(w,t) E*(w,t,p dw (111.2.11) 


Recalling the stationary case for which the mean square value 


is given ly 
.2 


E[y (t)] = / S (w) m 


( 111 , 2 . 12 ) 


— oo 


and comparing Eq.III-2.12 with 111,2.11 we can write 

Sy(w,t) ® S^. (w) K(w,t) E*(w,t) 

2 

= S..(w)|K(w,t)J' 


(III. 2. 12a; 


which is the time dependent power spectial density. 

Equation 111,2.11 after pecforming the integration gives the 
mean square response as 

y 

-2^w^t 

. [1. _L— y- (1- ^ ^ oos 2 + <irir5 2 ^ z 

2^ w 

^ o (111,2.13) 

Iq, 1 11,2, 13 can be further simplified by neglecting the 
higher order terms of , as ^ < < 1, we get the final expression, 

2 


y(t)’ 


- '"2Cw t It s*<o) 

D- . ] — ^ — 


2<, 


3 

’n 


( 111 . 2 , 1 ':-) 
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where . 

S.*(0) is the uniform power spectral density of white 

noise, 

III, 3 Numerical Approach 

The transient response of the single degree of 
freedom S 3 ratan, can also be evaluated, by finding directly the 
ensemble averages of the sample response functions. It may be 
written as 



(t) * E[y(t) y(t)J 



y^PyCy) 


dy 


(III. 3.1) 


where p-y(y) is the probability density function of the random 
variable^ T and is read as, probability of Y lying in the 
range y, y+dy. 

Discretization of the Eq. 111,3.1 gives, 

2 K 

<^y(t) =_S V^i^ (III.3.2) 

As the probability density function is not known 
before hand,c>~- (t) can be calculated by taking the arithmetic 

2 y 

mean of y^ as follows 

c — (t) * E — — — (111,3.5) 

• . : y i=1; : ' 

where i is taken across the ensemble. 
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Ho 111 in ear Response Ana^Ysis 

4 single degree of freedom systan of IXiffing’s. type is 
considered for the nonlinear response analysis. Ihe characteri- 
stic equation for such a system may be writ tea as* 

2 2 3 

y(t) + 2^y^y(t) + % y(t) + e w^y (t) = F(t) 

(III. 4.1) 

Ibr small 6 and for the white noise excitation theoretical 

that 

results are available and it can be shown the mean square value of 

y(t) is given by 
2 

^ ^ = 2-._ (III. 4. 2) 

y 

2 , , 
where <y~ is the mean square value of the response y(t) of 

o 

corresponding linear system. 

In the present thesis work the Monte-C5arlo technique is 
employed to investigate the effect of varying 6 on the system 
response characteristics, 

III, 5. Monte-Qarlo Response Analysis of an 
Aircraft Landing Gear Assembly to 
l axiing Induced Yibrations 

Finally, to demonstrate the usefulness of the Monte- 
Cjarlo approach, an idealised no del of the aircraft and its 
laaiding gear assonbly is considered and its response is studied, 
to the taxiing induced vibrations due to runway . rou^iness , 

Fig, III, 2 shows a sketch of the model. In this model, the 
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landing gear n^tss is neglected as it is very smll when 

compared to the mass of the aircraft, ibllowing are the equations 
of motion of such a system, ^or the derivation fief. [2^ may he 
referred. 

-k^(x^ - Xg) + ^2 ~ c^(x^ - x^) = y(t) 

(111.5.1(h)) 

where is the mass of the aircraft, c^ the damping coefficient, 
k^ is the stiffness of the landing gear, k^ the tyre stiffness, 
^1’ ^2 parameters of an equivalent linear system 

of the aircraft landing gear assembly. Bsr the details of the 
equivalent linearization technique employed, fief.jj^ may he 
referred. y(t) is the runway roughness the power spectral density 
of which is given by lq.II.4. 2. “ 

Eqs, 111,5*1 can be split up into three first order 
differential equations of the following form; 

^1 * ^2 

y, = ^ Ly(o - yJ (111.5.2) 

1 . 

y^ = 1/c., [kg y(t) - k^(y.j-*y 3 ) -Jc^y^ + C^y^ 
where y.j = y 'S^^- J and y^ = Xg (III. 5. 3) 

Init ial Qonditions At time t * 0, the syst^ is at 
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7 ^ = 0 

= 0 at t = 0 (III. 5.4) 



Once tile time history of y(t) is sinailated, the Eq.n,Ill,5*5 
can be integrated with the initial conditions given in 
Eq.*III,.5.4 and various response parameters such as, the f,m.$ 
values of displacement of the aircraft mass, the deflection of 
tyre, the velocity of the aircraft mass and the acceleration 
of the aircraft hbss can be computed. 

III.5 Results and Discussion 

The fallowing structural systems are considered:- 
1) 4 single degree of freedom system excited by band limited 

white noise. 2) 4 nonlinear system of Duffing*s type excited 
by band limited white noise 3) 4n idealized nKJdel of am aircraft 
and its landing gear subjected to taxiing and take-off induced 
vibrations due to runway rou^ness. The results are divided into 
two broad sections.!' In the first one the results of Monte-Qarlo 
approach are discussed and in the second one the results of 
transient response calculations based on ensemble averagesare 
presented. 

1) Monte Garlo Response Ahalvs is 

Two simple oscillators with critical damping ratios 
<^s= 0,1 and^» 0,01 are considered for the lin^r case, 3ample 
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tine history of the bard limited idxite input process I 0 
simulated digitally on the computer and response is computed 
by integrating the equation of notion using a fourth order 
ftmge^tttta integration sch«ne* The step l«»gtti of the integ- 
ration is chosen as one fourtieth of the undanqaed natural 
period of the oscillator. Integration is carried out for 
250 cycles of the response* Mean and Mean square values of 
the response are computed by taking the teiJporal averages 
over every 2*5 cycles* Results are presented in the form of 
a plot between the Mean square value and time in Pig,III*5, 

III, 4 * The dbtted line in the figures correspond to the 
stationary value of the response obtained from -Bie theoretical 
result. There is an exc^lent agreement between the similated 
and theoritical stationary values for casef* O.l. The inaximTam 
error is about 5*/* . Sox the Case*^* 0,01, there is a large 
differeiwe between the simulation and theoritical r^ults 
and the response takes significantly long time to attain a 
stationary value* 

Results of the nonlinear response analysis are presen- 
ted in Table III*1, Rig* XU. 5 shows the mean square r^ponse 

of the Duffing*s syst«Bs for foiir values of 6, - Also 

preset cd in the same figure is the result of a linear systaa 
for ccmparisonf The dotted line in each case represents the 
stationary value of mean sqtsare response obtained. by Bl.III.4.2. 
It nmy be infeirred from the result that as the amount of non- 
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linearity increases the mean square response reduces, and 
there is a sigaif leant reduction in the time taken by the 
response to reach stationarity, lor 6 * 0*2, the response 
has reached a stationary value after 15 cycles of response, 
whereas for € = 0*1, the system has not achieved stationarity 
even after 50 cycles of response. It may be concluded that 
the greater the amount; of nonlinearity the lesser will be the 
number of cycles of response needed for computing the temporal 
avearages » 


fhe mean square response in the trai^ient stage 
presented in the same figure does not isean anything as the 
transient response is non~stationary and an ensoahle average 
should be considered in place of temporal average. 


2. Transient Response Analysis 

IXiring the transient stage the response is not ergodic 
and an ensmble average should be considered in place of 
t^poral average. An ensemble of 36 sample response functions 
for the casef = 0,1 and 56 samples for the case 5= 0*01 are 
considered. Results are presented in Pigs. Ill, 6 and III .7. 


Also presented in the same figures are the results of 
Ghakro vort y^^ ^ for cou^ri so n purpo s es . The siimalatio n 


results are smooihed by taking neving average over a iCTfifch 
of 1 sec. or 5 cycles of response and plotted to the same 
scale. It is ci^r from the figure III.7 that there is a gDod 
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agreement between tbe tbeoritical and simulation result© 
for the case ^ = 0,1. However, for the case^= 0.01, the 
similation results are not so close to the theoritical results. 
The maximum a^or is about 25/* • For the case^ * 0.1 the 

response attains a stationary value very quickly, approxi- 
mat^y after 4 cycles. For the case^= 0.01, Iheoritically, 
the mean square r^ponse riches a stationary value after 
55 cycles, whereas the simlation results do not show any 
©tationarity in the range considered. To increase the accuracy 
more number of sample response functions should be considered 
and it is advisable to consider more number of cosine terms 
for the generation of sample input function. 



GH^m 17 

PIRS7 PASSi&E PROmSILIiTY AID OTHER SgATISTIOS 


17.1 Eor structural syistaiis operating in random environment, 
it becomes necessary to k 2 H>w certain statistics, such, as the 
first excursion time of a safe limit, 2) the amount of time 
spent above the safe domain and 3) the number of threshold 
crossings during the period of operation. Such statistics help 
in establishing the reliability and fatigue life of the 
structure. Bfeny authors jj3,14,15,l6,1'fJ have studied these 
problems and obtained solutions both analytically and numerically 
under certain basic assumptions. This chapter provides a dis~ 
cussion of the numerical approach to the pr obi on using the results 
of simulation outlined in Chapter II. The results are compared 
with those of the former authors. 

17.2. First Passage £robabilltv density 

First passage probability density, P^(T) may be defined 
as the probability that the random process -x(t) exceeds a 
threshold for the first time, during the time interval T to 
T+dt. 

A simple linear oscillator, excited by a white random 
process is considered in the following discussion. For such a 
system it has been proved faj assimiing poison process for liie 
arrivals of peaks, that 
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where a is expected number of crossings of barrier b and is 



(IV. 2.2) 


5 " is steady state root mean sqT:e.re value of the response. 
A is equal to 1 for b > > c' . 


is number of *0^ CTOssings. 

Ghandiramani and Go ok Litl Imve found out p (I) numerically 
from the statistics of safe intervals, shown in 1‘ig. IV.I, 

The excitation is simtfLated following the approach of Cook f2^ . 
The probability density function P^-C^) nay be written as 


) = P 


i .f • 


6 ( 1 ) + (1 - 




(IV. 2.3) 


where P. ^ stands for probability of instantaneous failure 
under the slationary starting conditions. Por *0* start condition 
clearly = 0 and Eq. IV. 2.3 may be written as 

P^(3:) = / 

In the present thesis work a tD* type barrier. }_I)ouble 
barrier^ 10 with *0* start condition is assumed. A sin^e degree 
of freedom linear oscillator subjected to band limited white 
noise is considered, A computer programme is developed whidb. 
cong)utes the response of the oscillator till it reaches, for the 
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first time, a threshold. lEae first i)assage time is stored in 
a dimensioned variable and the whole procedure is restaited, 
k separate subroutine simulates the random force. 1795 sucb first 
passage times are collected using five different random nianber 
generators. 

The probability density function is not unique as it 
depends upon the choice of end points of the interval for the 
range of the variable and the choice of the number E of class 
intervals in which the range is divided. The optimum values 
of E for different sample sizes N are given in Table 4.1 of 
Section 4.6. of ref. [2^ . E is chosen accordin^y. The proba- 
bility density function is estimted at the mid points of the 
E class intervals using the method outlined in Section 7.2.2, of 
the same reference. 

IY.5 Expected Amount of Time Spent Above a Threshold 

Consider a narrow band random process x(t), Eig. IIT.I, 

such as the response of a lightly damped linear oscillator to 

# 

white noise. The shaded area in the Eig. IV. 1 corresponds to 
the time spenb above a level ‘b*. It is r^uired to find out 
the ratio ^ ^ ~ which directly gives the expected amount 
of time spent per unit time above level *b* . This can be 
evaluated by monitoring the output of the single degree of ^ 
freedom system to a broad band excitation with the help of a 
digital computer. ‘ 
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Aml3rfcically it can be calculated as follows* Define 
a random process z(t) as 

z(t) s= 0 for x(t) < b^. 

- 1 foi\ x(t) ^ b 


(IV.3.1) 


How E Atj^ is given by 

t 

S At. = y = / z(t) dt (IF.3.2) 

^ o 

clearly y is a random variable for which the expected value 
can be calculated by 

t t 

E[y] = E / z(t) dt = / E{z(t)2 dt (IV’,3.3) 

o o 

1} ? 

= / ^ (0) X P[x(t) < ^] + 1 X p[x(t) > b] J dt 

“ / / PYCx,t) dx dt (IV *3 *4) 

0 0 


4s the system response is a Gaussian random process its pro- 
bability density function is given by 

• 2 ■ 

p^(x,t) s= — j:— V (-X / 2 o-^ (t)) 

V 27iS- W 


p 

where (t), transient mean square response and 

■ I3C 

m ^ v-g te at 

O O ?2n; G ^It} 


t 

/ enfc (z) dt 


^ 1 


(IV.3.6) 
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where 


'z ~ J 'b/2 5 " 


.(t) 


(I¥.3 .7) 


Eq.. IV»3»6 can be integrated using the standard procedures 
available for the numerical infcega?ation. 


I'7,4 Expected rate of crossings of a threshold : 

The expected rate of upcrossings of a level *b’ in a 
narrow band Gausian random process, such as the response of a 
simple linear oscillator excited by abroad band random process 
is given by 

oo 

v'*’ (b,t) =/*7f^ (h,^i;t) d^^ 


(17.4.1) 


Inhere f^^ (bj'^jt) is the joint probability density function of 
x{t) and xC'b) of the response. For the transient state it is 
given by 




2^ ‘~x ^x 


exp - 


_j 

2(1“ f^) '^xJTx 


(IT.4.2) 


where = b[x (t)J 

= e[x (t)3 


and f = BEx('t) ^(^)]/('^x "^x^ 
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Substituting Sq, 17.4.2 in Eq. 17.4.1 and carrying out the 
integration, it can be shown that 


/(t,t) = - 


Lxp( g(l+erfg) 

’ I 2 ) 


exp(- 


-j? 


2 - 


where g = 




Vx ^2C1- f''^) 


(iy.4.2) 

(IT.4.3) 


2 ^ ^o 

c- (t) = -f-T 

X 


. exp(-2^ I 


cos 2w^t 


1 -q 


t]j 


(iy.4.4) 


2 xS exp(-2'>ot) 

i (■‘5 = aTi" 


/ O 


1 


.2 


^ sin 2w^tj 

li |l - f^cos 2w^t -(-|l- Y 

sin 2»^t J ] 

(17.4.5) 


f(t) - 


K $ 


W 


O 






2hW 


-6^2 


exp(-2(; w^t) (l-cos 2w^t) J 


S. is the power 

d 

the 


spectral density of the input process and'y#^ 


mtural fr^uency of the oscillator. 


WnTno-rica l aoproa.^: 

Numerically it can be calculated as follows: 
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A sample record of consideiable iQigth is computed 
and is divided into a number of intervals. In each, interval 
the actual number of crossings of a specified lev^ »a* 
is counted and its ratio to the interval length gives 
directly the number of crossings per unit time. This will 
be a random variable as we have considered only one record* 

To get the expected value a number of such sample ^ould 
be considered and an aveiage across Ihe oisemble should be 
taken in a similar way in which the transient response is 
computed. In the present thesis work only six sample are 
considered* 

IV. 5 . Results and Discussion 

The response of a simple linear oscillator with 
critical damping ratio A = 0.1 is considered for the analysis 
of first passage and the related statistics. For the cal- 
culation of first passage probability 1795 sample response 
functions are considered. She safe domain is of double 
barrier type. A computer programme is developed, which 
automatically collects the first passage times, by noting 
the time at vhich for the first time the absolute value of 
the response exceeds level. Cta.ce it crosses the procedure 
is restarted and a new sample is generated. The total computer 
time taken for collecting 1795 samples is about two and half 
hours. The toaporal range is taken as 4 natural periods and 
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is divided into a large number of class Intervals, First 
passage probability density is calculated by deviding the 
number of samples falling in a particular interTal by the 
interval length and is plotted against the nondimensional 
time in Fig. 17.2. Results obtained by Crandall, (Siandi- 

ramani and Cook Dfl 

in the same figure, loland Yang also have obtained simlar 
results (3^3 • 

The results show a good agreement with the results 
of the previous authors except In the range. 4 2H to .64H^ 
where the probability densities obtained in this work differ 
significantly from the previous results. 

In Figs. IV.3. and IT.4 expected amount of tine 
.pent ahove 1 . and 2.- levels respectively by the oscillator 
response is plotted a^inst ncnd^ensionai time. The dotted 
line in the figures correspond to the theoretical result 

Tt Tir T-b is clear from the curves that 

obtained by the Eq. IV.3^S. It is c±ear 

.. j-b. no a-n excellent agreement between t e 

for1<r- level, there IS an exceixeii 

theoretical and simulation results approximately after 
,0 periods of response o»ards. during Ibe ttansi^t st^e 

simulation results predict sli^lF lesser values 

Tt may he said tl»t the expected amount 
theoretical values. It may he sa 

per unit time takes more time to attarn a 

of time spent 

trai ue for the case of sii»lation resu s# 

a+.ationary value lor 
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Itor 2 «-level, the disorepanoy between the theoretical 
and Bi^latlon reeuXts in «re pronounced. Sirnlatlnn reeulte 
sligl^tly overestingite the stationary value "tite nasi 
discrepancy being about 6-/. . However, in the tiansient 
region there is a large disorepanoy in the simulated and 
theoretical values, the simulated values being quite low 


compared to the theoretical results. 

In the H-g. IT.5, tiie erpeoted rate of uporossings 
of 1 ^-and 2 .- levels 1 in a narrow band random process 
The results show a similar b*aviour. 

The discrepancy in the simulated results is mainly 


etue to the inadequate sample sise. Only sir samples are 
considered for the calculations, theoretically the averages 

+vno PTLc^emble for infinite sample 
should be taken across the ensemOi 

fonotions, wherever the principle of ergodioity 

^Ild. this explains ol^rly the discrepancy between the 

j. ,w.n,otiT +<:5 in tile transi^t region, 
theoretical and simulated results m the 

where the response is not ergidic. 



CHAPTER V 
RAHim WAT.TT MOPEL 


V, 1 In the preceding chapters we discussed an approach to 
the stochastic process analysis in which the stochastic 
process of interest is simulated following its spectral rep- 
resentation and response characteristics are computed using 
the simulated time history. Throughout the analysis only the 
first and second moment functions are cGDg>uted. In the present 
chapter a different approach, the Markov process model of the 
stochastic process, is discussed. A brief description of Marlsov 
process approach is given below. 


¥.2 The Markov process is characterized by a distinct property, 
that given precise knowledge of the present, the future is 
independmt of past. Though this condition makes it highly 
restrictive, it happens that ngtny physical processes ei^ibit 
the Markov property. The response of a lin^r single degree of 
freedom oscillator is one such example, Doob C25j : first steted 
that the response, velocity, and displacement of 1iie lin^r 
system excited by white noise form components of Markov 
vector, Gau^ey has discussed its applicability to non- 
linear syst ems , Ghandiramani and Cook solved several 

first passage probloas u^ng the solution of Wang and Uhlcsnbech 
]j§] to the Pokker-Planck equation which dOiaract erizes the 
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Markove Process, Recently Poland |j^ proposed a random walk 
model which is applicable to a - sin^e degree of free- 
dom system. For a non-^hite excitation the solution is yet 
unsolved. Herein we describe briefly the Random Walk Model of 
!Eoland and extend it to a two degree of freedom systea 

subjected to white noise-excitation, 

IV. 3 I^ndom Wulk Model 


A general single degree of freedom systaa can be repre- 
sented as 

y + H(y>y) » J'(v) j y(o) *= o, y(o) * o (v,3.i) 


where P(t) is a sta-tionary Gaussian white noise with mean 
and correlation 


ElP(t3 * 0 

ElP(t^) PCtgl] » 2u 6(t2-t^) 


(V.5.2) 


6 is Dirac delta functioh. 

It is required to find 12ic- transioit conditional pro- 
bability of response y, y, at any timetv given that the system 
started from rest* 

This can be obtained by discretizing the Chapman- Kb lomogorov- 
Smoludiowski equation, for the two dimensional Markov process 
and may be written as 
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s[yi. yy V®-®-®] = s z r[v y„, 

X i) 21?!. yj . VV yn> Vi] 

(7.3.3) 

is the probability that the system is at 

m, 

^1* phase plane, at ..time tj^ given that the system 

is at rest to start with.^Tt^ is the second order transitional 
probability* 

B<i, 3«3 in^Jlles that all the possible states y^jj, y^ in -the 
phase plane, s-i have the possibility of r aching y^,. 

y. at time t, . So reduce the courplexity of the pro hi ©a, it is 
assumed tlmt once the velocity takes a definite value y^ than 
the displacement follows 111 e deterministic relation y. == y ^ t. .. 
The transitional probability function now reduces to 

The second assumption is that the velocities can chaaage only 
by small amounts. Consequently it is assumed that the velocity 
can either increase by one unit or decr^se by one unit, with 
probabilities p, and q. Now Bl. V,3.5 reduces to 
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Pi&i-V tj^_^/o,o,o]3: 

^C^l+V - ^1+1^*' Vlf'^'O^'3 (T.3.4) 

BqL*V,3.4 means that the structure which is at A presently as 
shown in Fig.T.1 was either at B and went up with probability 
p or it was at G and dropped down to A with probability 
The relations for p and q are given by 

p = ^ |j - H(y^, y^) At/^, (v.3,5) 

q * 1-p 

and ^ \ 271 Sq 6 t 

Eqs. V.3.4, and 7,3.5 characterizes completely the random 
walk model of the simple degree of freedom syst-em. However, 
for the computational conv lienee Eq. 7,3«4 iiB,y be written as 

P(m,r,s) = P(m-1, r-mf1,s-1) p(m-1, r-mfl, s-1) + 

+ P{W1, r-m-1, s-t] q(iin-1, r-m~1, s-1) (7,3«6) 

lAiej^e r, m, s are integer values representing discrete points 
in the phase plane as m - yj^r iiy - y^» s Zi t = tj^ 
and Ziy is defined as^*At* 
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Y,4. :First Tassage Probabili ty Density 

Random walk model can be applied to calculate tlie 
Plrst passage protatillty density, let us define p(y,y,t/0,0,0,r) 
or the prohahility that the response is '(y.y) at time t given 
that it is (0,0) at tbe start and all the probability masses 
are inside within the safe domain f. It can be shown that the 
probability density function of first passage time is giTen 

by p] 

g(t,r )=-'|t p(y»yiV 0 » 0 , 0 ,r) dy, dy 

^ (V.4.1) 

Numerically it can be calculated by the finite difference 

appr oxinei.t io n, 

PCt) - P(t ^ (V.4.2) 


g(t,r ) * C 


at 


where 

P(t + At) is the total probahility -mss at time t + M and 
P(t) is the total probability mass at time t. within the domain 

T.5. waiv Model for a ^o Degree of freedom Sygtg a. 

-p -FT.^^>,ion system may be represented 
A general two-degree of freedom sys-ceiu 

in lig ■^.3. can be represented by the following 
as shown in lag*'' *2 • J-*' ^ ^ 

equations of motion. ® ^2^ ** 

+ 02 (4-1) t = ?2(t) 

(v.p.D 
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where F^Ct) and independent white naise excitations 

with ^0* mean, Eq.s» V,6*1 are second order coupled simultaneous 
differential equations. It is assumed here that the equations 
can be uncoupled, using the norB&l mode approach Ihe un- 

coupled equations may be written as, 

yi + H(y^,y^) = f^(t) 

(y.5.2) 

7*2 + H(y2,y2) = £ 2 ^^^ 

• 2 

where H(y^,y^) - + w^ y^ 

♦ « 2 

H(y2, J 2 ) = 2 /■, 2 W2 J 2 + ^2 ^2 

f^(t) and f 2 (t) are white noise excitations having 
the following characteristics. 

S[f^(t)l = 0 , E[f2(tfl = 0 

E|i^(t), f^(t+73 * 271 6(7 ) (y.5*5) 

Ejf2(t) f2(t+T )] = 2x Sg 6 (t: ) 

It is required to find the transient probability of 
response y^» ^2 ^■hy time t, given that the systoa 

started from rest. 2^scretisation of the :• Chapman, Eomogorov - 
Smoluchowski equation for the above systen gives directly the 
transient probability of response as 
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p[S-ii. y^j. yzL- - 


Z £ Z S ^&1n> yio’ ^2p’ ^21’ 


n 0 p ^ 


/li^yil’^H’ ^Zk’ ^21’ V^lB’ ^1o’ ^21; ^24’ 


L4^yn>'’l3’ -Zk 


(T.5.4) 


Under the same assun5>tions as that of the single degree of 
freedom system, the fourth order transient probability function 
n can be reduced to a second order transient probability 


funotionj'i.g “'1 ‘’® written as 


■^4 ^^11’ ^13’ ^21t’ V^ln’ ^10’ ^2?’ =^21’ 

=J1, (y,„ y,^. v^n> yzp’ Vi> «Cyi3-yio-yih^*^ 


^ 1^21 ” ^2<l* 


...(■7.5.5) 


substituting Eq. 7.5.5 in equation 7.5.4 

^T* U T y ^ /0,0,0,0[ = 2 2 P 7^3 “ 7111^ 

^l^li’ ^13* ^2k’ ^21* nf / V ^ np 


• V -V Tit, t, \,/0, 0, 0, 0, 0[ xj(2[y-(j^j J’gk* k 
7o^» 721 yap^ ’ k-1 
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J\2 is a second order transition probability density 
function and gives tHe probability that velocity dianges 

from y^j^» ^ 1 i’ ^ 2 k ^k* . 

Now we assume that the velocity changes are likely to 
he in the closest nei^bourhood during a snail interval of 
time t. Consequently it is assumed that, in one time interval 
the velocity can either increase or decrease by one unit i.e. 

can become ^1134-1’ ^ ^ 1 n -1 

y^j^ can become yaiH-IV ^ 2 p -1 

=>i+n^pther 4 combinations associated with 

So there are altogether com 
these changes* 


p be -me probability of Increase of 
4 be the probability of decrease of y., 

^ be the probability of increase of 
be the probability of increase of y2 

n tbe transient probability now reduc® 
v> 4 4’ 

= 6(y-,ir ^2p<-lV^ 

+ 6 (y-» 


y ^tj t^^, 

-y 2p 



Substituting V, 5 *7 in V,5*6. 


p[yii, y^j. ygk. ya- ~ 

>^1/11+1’ ^13’ ^a+1 ’ '‘i'^2 

+ ^13* ^2k:+1’ Jgi’ t^O, 0,0,0, 

+ rlyii_i’ yi3> y2k-i’ ^a’ 

- 5[yii.i- ^ 13 ’ %c-r =^a- "1^2 


To calculate the one-step probabilities p^i> 2 ,^i 
use the statistics of the velocity increments 

<iiy-,> = P'1 ^ 

*2 * ^ 

<( Ay^)S > Pt£^ l> 

< ^ 72 ^ ~ ^2 ^2 ^2^"^ 2^ 

* ^ * 2 

Relating the statietlos to the system prahlems and 


(7.5.8) 

q.2 

(7.5.9) 

6X ^ 
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by £l ] 

<4y^> = y^^.) At 

^2 (V.5.10) 

< ^ y^ > = 2 % A t 

• 2 

< Ay^ > = 2 % S 2 At 

Eqns, V, 5*9 Sind Y«5#10. proTide four equations for the four 
unknowns P^ 4 P 2 #q-j Q'Ud. q^ which are found to be 

p^ = 1/2 0 - H(y^, y^) 

qi « 1-P^ 

.Pa = 2 D - HCyg, yg) (y.5-11) 

^2 

In the above formulation it is assumed that y^ and y 2 
are independent and their joint probability is equal to the 
product of the individual probabilities. 

As a check to the v^ocity of the above model, the 
recurrence equation Eq, V*5l8 in the limit as At ■* 0 can 
be shown to approach the PoMcer—Planck differential eqi^tion. 



V, 6 lirst Passage Probabilit y Density 

•Ph-r the Two Deeree of l^^eedoia System 


The random walk model derived above can be used "fco 
find the first passage probability density. Suppose we define 
the safe domain as < a. me coordinates x^, x^ 

of the original system can be expressed in terms of the 

generalized coordiiates by the following raations 


where 



(V.6.1) 

t and 


D is a matrix, 

Ihe oonditioisl picbablllty of the geheraUsed coor^tee 

(Here the oohditiehal sigh is droi^ea for eoiwehlehce) ^ 

How pLxi.^i.i2.^5l PD-I.yi.y2*y2^ 


We have 
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i I 


) 


1 ^2* 

1. 


-1 


' w ^ rj 



-1 - 


Henc 


e Pli^, i^, xgl = l®l 


pC y * 1^1 ^ * 1 


wliicli is known 9 


(V,6.5) 


(V.6.4) 


Now for ary "tinie ■t» 


pU. <aT*i 2&' t' y ^2! 

^ x; x| X. < a (T.6.5) 


Fir^ passage probability density. 

]p(!E) - PCT^-^tjJ 
At 


. g(ti n ) *= + 


(¥.6,6) 


where P(T) - prohaMllty »se at time t = I. 

P(It ^t)- Total protatllity ^es at time t = T^ At 

eroloeed toide the eafe domain defined h. < a. 


Y 7 'MnmPT-lGal AspSOtS 

some of the numerioal aspects a«h as the choice of 

n .cor At A* etc. and the numerical diff lenities 

i>ropsr ^ j? 4--u 

t»rea due to the maooraw limitation of the 

that may "be eiK^ountered 

QGiisqed in tbis section, 
computer are discussea 
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The phase plane is boimded in the y direction because 
of the natural boundary conditions 0 < p < 1, idiich means 

iH(y^y)] < ^ (T,7.1) 

At 

For a sxii5>ljfe linear os cilia "fcox Etj* V* 5 * 1 reducss to 

2^WQy+w^y< (V.7.2»a) 

r.P • 2 271 So 

- 2 ?w^ y - w^y i (7.7.2 b) 

Eq.. V. 5.2 represent a pair of strai^t lines as shown in 
Pig, V, 2* dividing the phase plane into a sub domain,. It is 
clear from the _eq.uations that only At need be specified for 
the numerical study, A proper choice should be made such that 
it gives a sufficiently large subdomain for a reasonable comr- 
putational accuracy with a desired stability. The stability 
here means a finite mean and a finite variance for all time 
[3 '7], k single degree of freedom system with critical damping 
ratio = 0.08, *= 2 units, and w^ = 5 cps is considered 

here for the study. Absorbing boun^-rios are set at ^ 1^ 
distance in the y direction for the first passage probability 
density calculations. Any probability mass diffusing across 
these barriers is considered lost and will not take par u iu. 
further diffusion, Thou^, the y direction is bounded to 
reduce the computational effort two artificial boundaries of 
absorbing type are set at 1«5 ^own in Pig. '7.2, 

In choosing this, it should be kept in 
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probability mass diffusing accross tfeese barriers should 
be of several orders of mi.|iiitude lesher than that lost 
accross displacenent barrier. 

Table V, 1 presents how the atemory requir^eisfes would 
vary with the choice of At. The number of cells into 
which Ihe phase plane is to be divided increases enomously 
with the reduction of ^t and soon suz^asses the mesnory 
limit of any computer. For the present work At = 1/t25th 
of the natural period of oscillator is tried, however the 
process is found to be still unstable. Results dhow a large 
deviation from the theoretical results, A minimum At of 
l/50C)th of the natural period of the oscillator has been 
suggested by the authors of the ^ndraa walk k) del* ^ cor- 
din^y the number of cells into which the fdmse plane is 
to be divided would be 120,000 and •fe.e maaiory needed would 
be approximately Ibree times the total number of cells i,e, 

36 OK, As the capacity of the existing computer is only 28K 
the problem can not be handled directly. 

Use of tape units to store "Sie value is also tried. 

But the time involved in the winding and unwinding operations 
is found to be too prohibitive. To overcome this sort of 
difficulty, the aulhors [2(3 Rave suggested an approach in 
which a number of cells amounting to approximately 
are grouped togeldier avd in the y direction the grid 
structure is unchanged. Probability mass diffuses into each 
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cell in each iiifcerval of time. Thus an averaging process 
occurs for the y direction. OJhis is not tried in the present 
thesis work. 

The problem for the case of two degree of freedom 
system is formulated but even attempts to solve this problem 
could not be cade because of the unwieldy coiaputer memory 
storage requiremoits which was felt even for the <^e of 
single degree of freedom system where the counter memory 
required was of the order of 100 E. 

ijr ■ , 

. ? • 



GHkPTM VI 
CQNCSLIBIOKS 


Sample functions of a stationary, Russian one 
one dimensional, and one variate random jaw cess, with a 
prescribed power spectral density are sinnilated ^gitally 
on the computer. Ihe samples are found to i^tisfy the t^ts 
for static narity. Also the power spectral density estimates 
from the saE5)les are found to agree with the target values. 
Parzen*s lag window is foimd to give a smoother and better 
estimate than Haniii4tg*s lag window. 

IChe transient and stationary response of simple linear 
oscillator oscillators with the dampii^ ratios ^ * 0.1 and 
S = 0,01 are computed. A moving average over a len|^of 
5 cycles of response is taken to smoolhen the. transient 
response. It is found tl&t for the case o'‘ — 0.1, the 

moving average values show a good agreement with the 
theoretically obtained results, Por the case ^ - 0,01, a 
larger sample size (only 50 ^mples areyconsidered in the 
present work) is recommended to reduce the discrepancy 
between the results obtained in the present work and the 
theoretical results. For the stationary response a sample 
record of i50 cycles of r^ponse is fo^d to be sufficient, 

L few nonlinear systeias of Duffing *s type are alro . 
considered for the response analysis, The stationary r.sus. 
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value of the response is coirqputed based on the temporal 
average over a length of 50 cycles* It is found tl^t with 
the increase of non-linearity there is a significant reduc- 
tion in the time taken to reach static narity* 

First passage probability density is found out 
numerically and is smoothened by taking a moving average 
over a length of 0*1 cycle. !Ihe results show a reasonably 
good agreement with those of previous authors. The expected 
amount of time spent above 1 tr'and 2^ level and the expected 

rate of upcrossings of lor'and 2 C* level, by the response of ' 
a simple oscillator are computed both analytically and 
numerically. The static mry values show an excellent 
agreement, 

A random-walk model for a two degree of freedom system 
is foimiula ted* However, it is found that Ihe meaory of the 
existing computer system is not sufficient to hasadle the 
problem. 
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U.BUE NO, 11,1 

RESDiaS OF TESTS POR ST4TIOMR3ZPY 


Level of Significance 

# 

o 

II 



S.No, Sequence 

Observed Yalues 

Theoretical bounds 


No. of 

No. of 

Ho. of 

No. of 


runs 

reverse 

runs 

reverse 



arrange- 

ments 


arrange- 

ments 

1 . Means 

22 

417 

15-26 

519-460 

2. Mean square 

values 

25 

586 




Note: The theoretical hoimds for l^e «mh«r of nins and 
the nunitoer of reverse arraDgemefits are taken from 
the tables 4-»11 and 4.12, pages 170-171, of 
ref erence |.22] • 
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Data Target power spectial density = 4 
Frequency limits 0 cps to 100 cps 
Maximum lag number m = 20. 


S.No. Frequency in 
cps 


£ 9 

Hanning’s 
lag Window 

parz^’s lag 
window 

4.13 

4.05 

4*06 

4.09 

4.06 

4.04 

4.16 

4.07 

4.18 

4.06 

4.13 

4*06 

4*05 

4.03 

4.09 

4.06 

4.24 

4.12 

4.41 

4.19 

4.52 

4.19 

4.13 

4.07 

3.83 

5*96 

3.88 

3.93 

3 .84 

3.91 


5.75 

5*59 

5^63 

3*72 

3*75 

3.78 


3.86 

3.82 

3.82 

3.85 

3.87 

3.88 



s 

Ifetuial frequency ^of tiie corresponding lin^r 
system is 5 .eps. = 0,1. 

Excitation power spectral density = 4 x 10'^ units. 


S.No, g Stationary mean square value of tlie response 




Simulation 

Theoretical 

1 

0.00 

10.50 

10.10 

2 

0.01 

7.54 

8.08 

3 

0,05 

5.77 

5,5 

4 

0.1 

4.41 

4.37 


5 


0.2 


2.97 


3.55 
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TkmE HD.V.1 


Systeci characteristics 

s= 0, 06, = 5 cps ~ ^ units 


S.No. At /f ^Av of 

0 j — j memrj loca- 

tions needed 


1 

,004 

50 

0.224 

5 

40 

2400 

2 

.002 

100 

0,1591 

8 

112 

10740 

3 

.001 

200 

0,1122 

10 

.317 

57100 

4 

.0004 

500 

0.071 

16 

1250 

360000 
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